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This  final  report  on  AFOSR  project  “Multilevel  Techniques  in  Large  Scale  Com¬ 
putation”  (AFOSR  number  F49620-92-J-0439)  is  described  according  to  the  following 
topics,  which  include  examples  of  how  these  applications  have  led  to  important  basic 
advances: 

•  First-Order  System  Least  Squares  (FOSLS) 

This  relatively  new  methodology  is  aimed  at  developing  a  well-posed  functional 
minimization  principle  for  a  broad  class  of  partial  differential  equations.  The 
goal  is  to  provide  a  general  formulation  that  is  elliptic  in  each  variable  so  that 
standard  discretization  and  multigrid  solution  method  yield  optimal  results. 
The  project  established  such  functionals  together  with  a  theoretical  foundation 
for  this  approach  for  diffusion-reaction  equations,  connection-dominated  flow, 
Stokes  or  Navier-Stokes  equations,  linear  elasticity,  Helmholtz  equations,  and 
problems  with  discontinuous  coefficients. 

•  Electromagnetics  and  Seismology 

New  multilevel  methods  for  modeling  very  general  wave  phenomena  have  been 
developed  that  involve  a  mixture  of  wave  and  ray  dynamics.  They  have  been 
shown  to  be  highly  efficient  on  realistic  electromagnetic  applications  involving 
homogeneous  and  inhomogeneous  media  for  various  irregular  exterior  domains 
geometries. 

•  Meteorology 

In  collaboration  with  J.R.  Bates  and  Yong  Li  of  NASA/Goddard,  we  developed 
a  multilevel  semi-implicit,  semi-Lagrangian  method  for  shallow  water  equations, 
based  on  potential  vorticity  and  divergence,  that  has  proved  to  be  more  efficient 
and  much  more  flexible  than  methods  currently  used  in  numerical  weather  pre¬ 
diction.  In  collaboration  with  Sam  Yee  of  the  Phillips  Laboratory  at  Hanscomb 
AFB,  we  have  developed  effective  multilevel  local  refinement  techniques  for  this 
approach. 

•  Multiscale  Atmospheric  Data  Assimilation 

The  main  bottleneck  today  in  atmospheric  and  oceanic  flow  calculations  is  the 
need  to  assimilate  into  them  a  continuously  incoming,  extremely  non-uniform 
stream  of  measurement  data:  optimal  assimilation  requires  the  solution  of  a 
huge  n  by  n  dense  system  of  equation.  Investigating  the  problem  over  the  last 
year,  we  have  found  that  multiscale  methods  can  solve  such  problems  in  just 
O(n)  operations.  Multiscale  data  assimilation  may  in  fact  be  even  faster  than 
direct  flow  simulation  itself. 

•  Molecular  Dynamics 

The  cost  of  simulating  macromolecular  dynamics  by  current  state-of-the-art 
methods  rises  very  steeply  with  problem  size.  Our  studies,  including  numerical 
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tests  on  simplified  models,  have  shown  that  this  steep  rise  in  cost  can  be  radically 
reduced  by  multiscale  algorithms. 

•  Robust  Three-Dimensional  Solvers 

In  collaboration  with  Uri  Shumlak  of  the  Phillips  Laboratory  at  Kirtland  AFB, 
we  have  developed  a  very  efficient  and  robust  three-dimensional  multigrid  solver 
that  is  well-suited  to  parallel  computation. 

•  Parallel  Adaptive  Refinement 

Work  has  continued  on  the  development  of  AFACx  for  more  general  applications 
(AFACx  is  a  multilevel  refinement  method  that  allows  for  effective  simultaneous 
computation  on  all  levels)  and  AMR++  for  a  broad  range  of  machines  and  ca¬ 
pabilities  (AMR+- 1-  is  a  C++  software  development  environment  that  supports 
general  code  production  but  provides  automatic  and  efficient  application  to  a 
host  of  advanced  computers). 

Before  the  mid  1980s,  multilevel  methods  were  slow  to  find  their  way  into  appli¬ 
cations  codes.  While  not  unusual  for  new  advances,  the  pace  for  multilevel  meth¬ 
ods  seemed  slower  than  one  might  expect.  One  of  the  likely  causes  was  the  begin¬ 
ning  misconceptions  about  multilevel  capabilities,  fed  by  frequent  failures  of  early 
developments.  However,  the  situation  today  is  drastically  different.  By  the  late 
1980s,  multigrid  methods  were  already  becoming  commonplace  in  codes  of  certain 
application  areas  and  they  were  rapidly  appearing  in  others.  For  recent  develop¬ 
ment,  see  the  information  and  codes  accessible  through  MGNet,  especially  papers 
from  the  Sixth  and  Seventh  Copper  Mountain  Conferences  on  Multigrid  Methods: 
ftp://na.cs.yale.edu/pub/mgnet/www/mgnet  or  http://na.cs.yale.edu/mgnet/www/ 
mgnet.html.  It  is  worth  noting  that  the  computational  approaches  developed  over  the 
past  several  years  in  our  research  on  high-performance  CFD  (introducing  advanced 
nonelliptic  multigrid  methods)  are  in  the  core  of  a  major  project  now  being  proposed 
by  several  senior  investigators  (headed  by  Dr.  Jerry  South)  at  NASA/Langley.  Their 
purpose  is  to  develop,  by  year  2000,  multigrid  solvers  for  3-D  aerodynamics  that  will 
be  two  orders  of  magnitude  faster  than  current  solvers. 

An  important  ingredient  in  the  effective  progress  of  an  evolving  discipline  is  the 
organized  dissemination  of  developing  knowledge.  The  project  team  has  been  actively 
involved  in  such  technology  transfer,  starting  with  the  Multigrid  Workshop  at  ICASE 
in  1978.  To  date  in  the  U.S.  alone,  the  proposers  have  organized  seven  Copper  Moun¬ 
tain  conferences,  three  major  workshops,  four  short  courses,  and  eight  tutorials.  In 
addition,  the  proposers  have  developed  texts,  guidebooks,  bibliographies,  videotaped 
lectures,  and  curriculum,  and  have  acted  as  a  resource  of  information  and  informal 
consultation. 

There  are  many  aspects  of  the  current  research  project  that  have  influenced — 
and  have  been  influenced  by — the  work  of  others  through  collaboration  and  informal 
discussion.  For  example,  there  continue  to  be  many  short  term  visitors,  institutional 
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affiliates,  and  informal  collaborating  scientists  who  have  interacted  with  research  and 
development  aspects  of  the  project.  Moreover,  there  were  two  post  doctorates  and 
two  sabbatical  visitors  who  worked  directly  on  project  activities.  While  such  support 
was  not  a  formal  part  of  the  AFOSR  funded  project,  this  mutual  interaction  with 
collaborators  was  essential  to  its  relevance  and  progress. 

Further  Information 

More  information  on  these  topics  (including  algorithm  details,  theoretical  and  nu¬ 
merical  results,  practical  applications,  and  historical  and  philosophical  perspectives) 
can  be  obtained  via  anonymous  ftp  at 

amath.colorado.edu/pub/multigrid 


and 

amath.  Colorado  .edu  /  pub  /fosls 
or  via  the  World  Wide  Web  at 

http://amath-www.colorado.edu/appm/faculty/stevem/Home.html. 
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